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Abstract. We develop a new mathematical model for describing a dynamical 
system at limited resolution (or finite scale) , and we give precise meaning to the 
notion of a dynamical system having some property at all resolutions coarser 
than a given number. Open covers are used to approximate the topology of the 
phase space in a finite way, and the dynamical system is represented by means 
of a combinatorial multivalued map. We formulate notions of transitivity and 
mixing in the finite resolution setting in a computable and consistent way. 
Moreover, we formulate equivalent conditions for these properties in terms of 
graphs, and provide effective algorithms for their verification. As an application 
we show that the Henon attractor is mixing at all resolutions coarser than 10~ 5 . 



1. Introduction and statement of results 

The theory of dynamical systems has experienced tremendous growth and de- 
velopment throughout the recent decades. Nevertheless, in spite of the undeniable 
importance of the results describing complicated dynamical properties many of 
these results are abstract existence or genericity results. This can be problematic 
for example in applications, where generally more concrete and quantitative results 
are desirable, especially if one is studying the properties of some very specific objects 
for which analytic methods fail. It is therefore natural to try to use numerical and 
computational techniques, and indeed there is an enormous amount of literature 
in this direction. Unfortunately, this approach is limited to finite resolution and 
to finite time and is therefore unsuitable for direct rigorous verification of infinite 
time or infinitesimal scale properties without additional theoretical interpretation 
of the results of computations. 

This combination of numerical methods and theoretical interpretation can indeed 
be remarkably successful in certain situations. In many cases non-trivial mathemat- 
ics can reduce the problem to the verification of certain "inclusion conditions" or 
the existence of certain geometric structures which can be rigorously verified numer- 
ically; then this information can be used to develop symbolic dynamics or prove ex- 
istence of certain trajectories, e.g., heteroclinic connections. This approach includes 
the application of some topological methods, like the fixed point index or the Conley 
index. We mention for example [Illl[Il[TMTg|15H23^ and 
refer the reader to references therein for further information. However, the bottom 
line is that by the intrinsic limitations of numerical computation, only "robust" 
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phenomena, i.e., phenomena which persist for small perturbations of the system, 
can be proved in this way. Indeed, the requirement that all estimates be rigorously 
bounded, necessarily implies that there is always a little margin of perturbation 
in which they continue to hold true; sec [30 for an interesting discussion of the 
properties which can be rigorously proved numerically, referred to as "inheritable 
properties" . This is of course in many ways also one of the strengths of the methods 
rather than just a limitation, but it does mean that "unstable" phenomena, i.e., 
dynamical features which exist for one system but perhaps not for nearby systems, 
are essentially undecidable, both in a theoretical as well as in a practical sense: 
see [J] for a formal discussion of this problem. The best we can hope for in this 
situation is to be able to prove that such phenomena occur for "some" , or per- 
haps even "many" , systems close to the one of interest, but even this strategy is 
sometimes extremely hard to implement. 

One specific example of this situation is the famous and very well studied Hcnon 
family of two-dimensional diffeomorphisms given by 



This family was introduced by Henon in 15! where he focused particularly on the 
parameter values a — 1.4, b = 0.3, now commonly referred to as the "classical" 
parameter values. For these parameter values, non-rigorous computer simulations 



Figure 1. A numerical approximation of the attractor of the 
Hcnon map for the classical parameter values. 

suggest that most orbits converge to a transitive attractor with a complicated frac- 
tal geometric structure (see Figure [T]), and exhibit chaotic dynamics. While the 
existence of an attractor follows by relatively elementary arguments, it turns out 
that it is extremely difficult, and indeed arguably impossible, to make any rigorous 
assertions about the infinitesimal structure of the attractor. Indeed, the geometry of 
the Henon map suggests the occurrence of tangencies between stable and unstable 
manifolds, and in fact in [3] it is proved that for b sufficiently close to 0.3 there ex- 
ists some parameter a € [1.392419807915,1.392419807931] for which a homoclinic 
tangency occurs (this is proved using some clever application of Conley index the- 
ory combined with some numerical "inclusion" estimates as discussed above) . Such 
tangencies are a well known source of instability and bifurcations and associated 
to a variety of dynamical phenomena such as infinitely many sinks [34 , 36 , "small" 
stochastic attractors [TU1I55] and others [2] . In particular the dynamics is extremely 
unstable and there seems to be no hope to establish the actual dynamics of any 
single given parameter. For small values of |6| it is nevertheless possible to show 
that there is a positive Lebesgue measure of parameters a for which the dynamics 
is stochastic in a very well defined sense, thus showing that stochastic dynamics 
has in some sense positive probability. This was first proved for b = in [18] . with 
some generalization to other one-dimensional maps in 123 [3 5, 44, 45], and for 
b 7^ in [SJ[5], with some generalizations in 33,47]. A quantitative approach was 



H a , b (x,y) = (l + y - ax 2 ,bx). 
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developed in |25j with some explicit measure bounds. Unfortunately, further ex- 
tensions of these results to larger values of to include for example the classical 
parameter values, seem for the moment still out of reach. We also emphasize that, 
as mentioned above, these results still only provide a probabilistic description of 
the parameter space, and do not in general yield information about any particular 
given parameter. 

Our goal in this paper is to approach the problem from a radically different 
point of view, which provides a different kind of information and is hopefully much 
more flexible and versatile, and more widely applicable. The underlying philosophy 
is that only finite precision is usually of practical meaning for applications, and 
all phenomena which take place below certain resolution are of no real importance, 
either because of the limited accuracy with which a mathematical model describes a 
physical system of interest (due to noise or truncations, for example), or because of 
the finite precision of measuring devices which provide observable properties of the 
system. Our goal is therefore to sketch the beginnings of a theory of finite resolution 
dynamics and to give an example of a non-trivial implementation of the concepts 
developed in such a theory. The starting point is defined by two basic "axioms" or 
"principles" which finite resolution dynamics should verify: 

(A) Computability and (B) Consistency 

These principles can in practice be implemented in a variety of ways, but they 
embody the fundamental ideas that: (A) dynamical properties should be formulated 
in such a way that they are rigorously verifiable using computational methods, in 
particular they should be related to finite time and finite resolution properties of 
the system; (B) computations at a finer scale should yield more valuable results in 
the sense that if a dynamical property holds at some finite scale then it should also 
hold for any coarser representation of the same dynamics. 

As an example of the application of these ideas we define a notion of combinato- 
rial mixing which, in the framework of finite resolution dynamics, is analogous to 
the standard notions of topological or measure-theoretical mixing for topological or 
measurable systems, respectively. We shall prove the following 

Theorem 1. The Henon attractor is mixing at all resolutions > 1CP 5 . 

We emphasize that, once the appropriate definitions have been put in place, this 
is a non-trivial rigorous statement about the dynamics of the Henon map for the 
classical parameter values, obtained by computer-assisted methods. 

The structure of the paper is as follows. In Section [5] we introduce the defini- 
tions which are at the foundations of our theory of finite resolution dynamics, and 
we state and prove our main abstract result. More specifically, we give a precise 
formulation of Axioms (A) and (B), define what we mean by the statement that 
a certain abstract property is verified at all resolutions coarser than some e, and 
demonstrate the key fact that this can be rigorously proved by computer-assisted 
arguments. In Section [3J we formulate natural notions of transitivity and mixing in 
the finite resolution setting and prove that these finite-resolution definitions satisfy 
the required consistency conditions (B). In SectionUwe provide explicit algorithms 
for the verification of transitivity and mixing at finite resolution, and thus in par- 
ticular provide a constructive proof that they satisfy the computability conditions 
(A). In Section [5] we describe the particular application of our theory to the Henon 
map and we prove that this map is mixing at some finite but very fine scale. 

This paper is accompanied by efficient and flexible software programmed in the 
C++ language and made freely available at [35] together with raw data for the 
applications discussed in Section [SJ 
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2. Abstract theory 

One possible approach to represent a dynamical system in a finite way is based 
on partitioning, that is, subdividing the phase space into a finite family of compact 
sets with nonempty and disjoint interiors. This approach goes back several decades 
and has proved extremely successful for describing certain classes of systems such 
as uniformly hyperbolic diffeomorphisms and flows [7HH1I1I] and has been recently 
further developed and used in various applications (see 32, 37, 42] for some ex- 
amples). However, working with this kind of discretization in more general systems 
may not well reflect the actual topology of the phase space. Therefore, in this paper 
we introduce the idea of using a structure more closely related to the topology of 
the phase space: an open cover. As we shall see below, the non-trivial overlap be- 
tween elements of an open cover provides a crucial ingredient for the development 
of a formal theory of finite resolution dynamics. 

2.1. Open covers. Here and for the rest of the paper, X will always denote a met- 
ric space; without loss of generality, we shall assume that the metric is bounded. This 
generality is convenient, because in applications, like those considered in Section[5l 
it may be necessary to work with the space X being a bounded, not necessarily 
open, subset of M" with the induced metric. An open ball centered at a point x £ X 
and of radius r > will be denoted by B(x,r). For a subset U C X, its diameter 
diam U is the supremum of the distances between any two of its elements. Since X 
is bounded, the diameter of any subset of X is finite. For a family A of subsets of 
X we denote their union ^y \A\. 

Definition 2.1. A finite family U of open subsets of X such that X = \U\ is called 
a cover of X. 

We are going to use the elements of U as a finite approximation of the topology 
on X. In general, the cover IA provides a better approximation if it consists of 
smaller elements. 

Definition 2.2. The outer resolution 

K + (U) := max{diam(?7) : U £ U] 
of a cover 14 is the maximal diameter of its elements. 

Intuitively, if two points x, y £ X are at a distance larger than the outer resolu- 
tion of a cover then they are well distinguished from each other by the cover since 
they must belong to distinct elements of the cover. 

Remark 2.3. For a totally bounded metric space (e.g., a bounded subset of the 
Euclidean space with the inherited metric), there exist covers with arbitrarily small 
outer resolution. 

Definition 2.4. We say that a cover U\ of X is finer than another cover U2 of X 
(or, equivalently, IA2 is coarser than IA\) if every element of U\ is contained in some 
element of IA2, and every element of U2 contains some element of U\. 

Note that in this definition we require a mutual relation between the elements of 
U\ and U2 in order to ensure that IA2 is indeed a coarser description of the topology 
of X and does not recognize any finer topological structures of the space than U\ 
does. This relation between (essential) covers defines a partial order on the space 
of (essential) covers of X, so we can write U\ -< U2 if U\ is finer than U2- 

In this paper we shall focus on essential covers, which are in some sense minimal, 
that is, no element can be removed from them without leaving a considerable part 
of X outside the cover. 
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Definition 2.5. A cover U of X is essential if there exists e > such that every 
U £ U contains some point x £ X such that B(x, e) C U and B(x, e) H W = for 
all WeW\{?7}. 

Notice that if a cover is not essential, i.e., there is an element U of the cover that 
does not satisfy the condition above for any e > 0, then U must be contained in 
the union of closures of the other elements of the cover, and thus is in some sense 
superfluous; therefore, covers that are not essential may introduce a false feeling of 
the topology of the space, and indeed, as will be made clear in the sequel, some 
features of our theory do not go through without this assumption. Because of this 
reason, from now on we shall always work with essential covers. 

Remark 2.6. From any cover it is possible to create an essential one (without 
increasing 1Z + , but possibly decreasing the number of elements); see Appendix lAl 
for details. 

Remark 2.7. For essential covers, the definition of finer and coarser covers given 
in Definition 12.41 reduces to the following: A cover U\ of X is finer than another 
cover Ui of X (or, equivalently, U2 is coarser than U\) if every element of U\ is 
contained in some element of lAi. Indeed, this automatically implies the second 
part: every element of U2 contains some elements of U\. To see this, just suppose 
by contradiction that there exists some U2 £ U2 which does not fully contain any 
element of U\. Then, since every element of U\ is contained in some element oiU%, 
it follows that U2 \ {U2} is still a cover, contradicting the fact that it is essential. 
We thank Francesca Aicardi for this observation. 

2.2. Combinatorial maps. We shall use the symbol — o to denote a possibly mul- 
tivalued map between two sets. Let A be a bounded metric space and U an open 
cover of X. 

Definition 2.8. A combinatorial map is a multivalued map J-:U—°U. 

Since U is finite, J 7 is a finite object that can be represented in a purely combina- 
torial way. If / : X — > X is a map, we can use combinatorial maps to approximate 
the map /. 

Definition 2.9. We say that a combinatorial map J-:U^Uisa. representation 
of a map / : X — > X if for every U £ U we have 

F{U) D{W£U:Wn f(U) ^ 0}. 

The ideal situation would be to have a representation where we have an equality 
•F{U) — {W £ U : Wr\f(U) 7^ 0} but computing such a representation is generally 
not possible in practice due to numerical approximation and computer round-off 
errors in the calculation of a guaranteed outer bound for f(U): see Remark [2.141 
below. At the other extreme we have the trivial combinatorial map which maps 
each U to all elements of the cover. This is a combinatorial representation of / but 
gives no information whatsoever about /. To control in some way this issue, and 
keeping in mind also that generally J 7 gives a better approximation if it is defined 
on a finer cover, we introduce the following 

Definition 2.10. The outer resolution 

K + (T) := max{diam[/,diam|J"(L/)| : U £ U} 

of a combinatorial map J- : hi — o U is the larger of the maximum diameter of the 
elements of the cover U and of the maximum diameter of their images. 

If / is a continuous map with some Lipschitz constant L > and the elements 
of U are of similar size (e.g., balls of the same radius), then in principle one should 
expect TZ + {F) ~ max{l, L}1Z + (U) for a reasonable representation T of /. 
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Definition 2.11. Let U\ and U2 be two covers of X. We say that T\ : U\ — o U\ is 
finer than Ti'.Ui —° Ui (or, equivalently, T% is coarser than if Wi is finer than 
U2 and for every U\ £ U\ and every U2 € % such that C/i C C/2, every element of 
F\(U\) is contained in some element of Ti^U-i)- 

Intuitively, T\ is finer than T2 if the images of T\ are smaller than those of T%. 
Indeed, it follows immediately from the definition that if T\ is finer than T2 then 
\T\(U{)\ C 1^2(^2)1 for all Ui £ U x and U 2 £ U 2 such that Ui C U 2 - A key point of 
our approach is to develop some dynamical notions which are actually computable. 
In particular, to take advantage of the notions of combinatorial maps we need to 
be able to compute such maps. 

Definition 2.12. / is e- computable if there exists a cover U of X and a method 
for computing a combinatorial representation J-:U^>Uoif with 1Z + {J : ) < e. 

This essentially depends on how explicitly we know the map / and is generally 
not a serious issue: see Remark 12.171 We conclude this section with a series of 
remarks concerning the definitions above. 

Remark 2.13. Unlike in the case of covers, the relation of being finer does not 
define a partial order between combinatorial maps. Consider the following example: 
X = {a, 6,c} with the discrete topology, U\ — {{a}, {b}, {c}}, U2 = {{a,6},{c}}, 
U 3 = {{a,b},{b,c}}, Ty. {a} ^ {{a}}, {6} {{a}},{c} k> {{ c }}, T 2 : {a, b} ^ 
{{a,b}},{c} -> {{c}}, J3: {a, b} 1 ^ {{a,b}},{b,c} ^ {{b,c}}. Obviously, U x < 
U2 -< U3. Moreover, T\ is finer than T<x and Ti is finer than ^3. However, T\ is not 
finer than J" 3 , because for U\ := {b} C J7 2 := {b, c} the element {a} £ Fi{U\) is 
not contained in any element of ^(L^) = {{b, c}}. 

Remark 2.14. Every map /: X — > X has a unique minimal representation relative 
to a given cover U defined by 

T f ,u(U) :={W£U:WC\ f(U) ^ 0}. 

This is the obvious "abstract" definition of a representation of / which would be 
natural if it did not have to be explicitly and rigorously computed. It is minimal 
in the sense that the image of any U £ U by any other representation T of / 
contains T^uilJ). We emphasize, however, that although some representation T 
of / can usually be computed in a relatively straightforward manner, the minimal 
representation for the same map need not be computable at all. This is because, in 
general, when using rigorous numerical methods, the best one can typically compute 
is some outer approximation b\j of the image f{U) of each U £lA. Therefore, it 
is usually possible to compute some WjJ containing all the elements of U which 
intersect By, good for constructing a representation of /. However, in general it 
may be impossible to determine which of the elements of W(j actually do intersect 
f(U) and which do not. This justifies our definition of a representation of / in which 
we do not require that it is actually the minimal representation, and we allow that 
■F(I7) contains some superfluous elements of U. 

Remark 2.15. The composition of combinatorial maps can be easily defined. For- 
mally, let 2 U denote the set of all possible subsets of U and define the map T: 2 U -> 
2 U as T(A) = UueA^iU) for any A £ 2 U . Then for all n > 1 we can write 
jrn _ jm ^ where i-u -^2 U is the "embedding" i(U) := {U}. In the same way 
we can define the composition of two different combinatorial maps J- and Q as long 
as they are both defined on the same cover. It is easy to see that if J- and Q are 
representations of /, g : X — > X, respectively, then Q o T is a representation of g o /. 
In particular, if J 7 is a representation of / then T n is a representation of f n for 
every n > 0. 
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Remark 2.16. Note that if T\,Q\:U\ — ° U\ are finer than Ti,Qi \ U2 —° Ui-, re- 
spectively, then Q\ o T\ is finer than Q2 o Ti- In particular, if T\ is finer than T<x 
then also the same relation holds for iterations of these maps, that is, J-"™ is finer 
than Indeed, take any U\ £ U\ and any U% £ U2 such that U\ C U%. Take any 
Wi £ Si(J"i(?7i)). There exists V\ £ 7i(J7i) such that Wi G C/i(Vi). Since .Fi is 
finer than J^, and U\ C E/21 for this V\ £ F\{U\) there exists some V2 £ ^2(^2) 
such that Vi c V2. Since £1 is finer than G2, and Vi C V2, for this W\ £ Qi{V\) 
there exists some W2 £ £2(^2) such that W\ C W^2- As a consequence, Wi C W2, 
where W% £ Q2 [F2 (^2))- Note that this proof does not go through if we replace 
"and every U2 £ U2" with "there exists C/ 2 £ U2" in Definition 12. Ill 

Remark 2.17. We make some brief remarks on the computability of combinatorial 
maps. If /: M" — > M" is a continuous map defined by means of an explicit formula 
involving only elementary operations (addition, multiplication, etc.) and simple 
arithmetic functions (like the trigonometric functions, <Jx, or e x ), and an open 
bounded set X C M. n can be found such that f(X) C X, then one can compute a 
finite representation of / : X — >• X using the concept of interval analysis [33] . For a 
cover U of X that consists of products of open intervals, a set B\j containing f(U) 
can be computed for each U £U, where B{j is also a product of open intervals. Then 
we define W[j to be the union of those elements of U which intersect B{j (such a 
set can be computed easily). The multivalued map f:M3[/i-> w(j C U is then a 
combinatorial representation of / on X . Obviously, the smaller the elements of the 
cover IA are taken, the smaller the outer resolution of the constructed representation 
should be expected. This argument proves that a finite representation of dynamics 
for a wide class of maps is computable at virtually any resolution. We shall use the 
idea described above in Section [5] for the analysis of the Henon map. 

2.3. The axioms of finite resolution dynamics. We are now ready to formalize 
the axioms of computability and consistency introduced in Section [TJ based on 
the notion of combinatorial maps. Let V :— ^{jF) be a predicate concerning a 
combinatorial map J- . 

Definition 2.18. We say that V is a finite resolution property if it satisfies the 
following two conditions: 

(A) Computability: V is computable, that is, there exists an algorithm 
which for any combinatorial map F can establish in finite time whether 
T satisfies V or not. 

(B) Consistency: V is consistent, that is, for any pair of combinatorial maps 
T\ and F2 such that T\ is finer than F2, we have V(J-i) V(J-2)- 

We remark that the definition is given purely in terms of combinatorial maps 
with no reference to the underlying space or any map whose representation J- might 
be. Thus we have an a priori notion of what it means for a property V to be an 
"acceptable" property for the investigation of finite resolution dynamics. In general 
we will apply this definition to combinatorial maps which arise as representations 
of some particular map / : X — > X on a metric space X and thus obtain some 
coherent statements about the "finite resolution" dynamics of the map /. 

We make here a simple observation which is crucial to motivate the theory. Part 
(B) in the definition above is not in itself sufficient to define what it means for the 
property V to hold at "all resolutions > e" . Indeed, assuming that V holds for some 
combinatorial representation J- with respect to some cover IA with outer resolution 
1Z + (J-) < e, it is easy to construct another combinatorial representation T with 
1Z + (JF) > e but which is not coarser than T and therefore is not automatically guar- 
anteed to satisfy V . In other words, choosing a cover with a larger resolution does 
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not necessarily guarantee that the cover is coarser, and therefore the consistency 
condition does not provide any information. We need some notion which guaran- 
tees that V holds for all combinatorial maps at sufficiently coarse scales, based on 
the fact that V has been verified at some scale for a single particular choice of 
combinatorial map. To formulate this notion we need to introduce an additional 
definition. 

2.4. Inner Resolution. The following definition and its application constitute in 
some sense the key idea of this paper. 

Definition 2.19. The inner resolution of a cover Li is 

n~{U) := sup{d >Q:VxeX3UeU: B{x,d) C U}. 

The inner resolution of a combinatorial map J- : Li — ° Li is 

K-{T) := TL-(U). 

Positive inner resolution is a crucial feature that distinguishes our approach from 
one based on partitions, where the inner resolution is 0. 

We have used the convention here that B(x, 0) = {x} so that the inner resolution 
is the supremum of the numbers d > such that every ball B(x, d) C X is contained 
in some U £ Li, or if such d > does not exist. In a connected space, the quantity 
1Z~ (Li) can be interpreted as the minimal width of overlapping between adjacent 
elements of U. Intuitively, if the distance between two points x,y £ X is smaller 
than the inner resolution of a cover then these points can be identified by means of 
belonging to a common element of the cover. 

Remark 2.20. If X is compact then we can apply Lebesgue's number lemma to 
any open cover U of X in order to know that there exists a number 8 — 8 (Li) > 
such that every subset of X whose diameter does not exceed 8 is contained in some 
element of the cover. Then obviously TZ~(Li) > 8/2 > 0. 

Remark 2.21. In certain cases, e.g., if X is a bounded subset of a Euclidean space, 
for any e > 8 > it is possible to construct a cover Li of X such that lZ + (Li) < 2e 
and 1Z~ (Li) > 8. In particular, this is true if there exists a finite number of points 
{xi, . . . , x n } C X such that Uq := {B(x i7 e — 6) : i = 1, . . . , n} is a cover of X. 
Indeed, let us consider Li := {B(xi,e) : i = 1, . . . , n}. Obviously, U is a cover of X 
and 1Z + (U) < 2e. Moreover, every x £ X is within the distance of r < e — 8 from 
Xk for some k £ {1, . . . , n}, and thus B(x, 8) C B(xk,e) by the triangle inequality, 
soK-(Li)>8. 

Remark 2.22. The inner resolution also provides a lower bound for how "thick" 
the cover elements are in terms of containing a ball of a big enough radius. Define 
T(W) := sup{d > : VZ7 £ Li 3x £ X : B(x,d) C U}. (See Figure 1 for an 
elementary example.) It is easy to see that for an essential cover Li of X we have 
1i-~{Li) < T(Li). Indeed, take any positive number d < lZ~(Li). Take any U £ Li. 
By the assumption that Li is essential, there exists x £ U such that for some e > 
the ball B(x,e) is contained in U and disjoint from all the other elements of the 
cover U. Since d < 1Z~(U), for this particular x there exists some U' £ Li such 
that B(x,d) C U' . However, U is the only element of Li that contains x, so clearly 
U' = U, and thus B(x,d) C U. This shows that for every U £ Li there exists an 
x £ Li such that B(x,d) C U. Since d < lZ~(Li) can be arbitrarily close to lZ~(Li), 
the supremum of these numbers is greater than or equal to lZ~(Li). 
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FIGURE 2. An illustration of the quantities TZ + (U), TZ~(U) in Def- 
initions and and T(U) in Remark l2.22l for a sample cover 
U = {U 1 ,U 2 ,Uz,U i }. 

2.5. Finite resolution properties for all resolutions > e. We are now ready to 
state precisely what we mean when we say that a property holds at all resolutions 
coarser than some e, and to state and prove our main abstract theorem. Let V be 
a finite resolution property and let / : X — > X be a map for which there exists a 
representation whose outer resolution is < e. 

Definition 2.23. / satisfies V at all resolutions > e if V is satisfied for every 
representation F:U^Uoif with > e. 

We are using here the notion of inner resolution to quantify the resolution of a 
cover, but recall from Remark l2.2f I that in many cases the inner and outer resolution 
can be chosen arbitrarily close (up to the factor of 2). Notice also that the definition 
requires V to hold for any cover IA and any combinatorial representation T : hi —° U 
as long as 1Z~(J-) > e. ft is therefore a priori an unverifiable condition. However, 
we have the following 

Theorem 2. Suppose V is a finite resolution property. If V holds for a represen- 
tation J-"o of f with TZ + (J 7 q) < e then V holds for all representations T of f with 
1Z~ (J 7 ) > e. Thus f satisfies V at all resolutions > e. 

Proof. We present the proof in the form of two lemmas. 

Lemma 2.24. Let U\ and IA2 be two covers of X . If lZ + (Ui) < 1Z~ (U2) then 
U X <U 2 . 

Proof. Take any number r such that 1Z + (Ui) < r < Tlrili-i)- Let U\ & U\. Consider 
a ball B(x,r) for some x € U\. From the fact that r > TZ + {U\) it follows that 
U\ C B(x,r). Since r < "RrifAi)-, there exists Us G U2 such that B{x,r) C U2, and 
thus Ui C U%. Now consider any U2 £ ^2- Then (by Remark I2.22p there exists a 
ball B(x,r) C Us- Since IA\ is a cover of X, there exists U\ £ IA\ such that x G U\. 
Since dianal^ < r, we have Ui C B(x, r), and thus Ui C U%. □ 

Lemma 2.25. IflA\ andlA2 are two covers of X , T\ : IA\ —° U\ is a representation of 
some map f : X — > X, and 1Z~ {IA2) > lZ + {J-\), then any representation T2 ■ U2 —° 
U2 of f is coarser than J-\. 

Proof. By Lemma 12.241 we have IA\ ~< IA2 ■ Take any U\ € U\ and any U2 € IA2 
such that U\ C U2- Take x € U\. Take a number r such that TZ + (J 7 i) < r < 
r Rr(JA2). Then there exists W2 G U2 such that B(f(x),r) C Y\?2- Since T2 is a 
representation of / and x E U2, ^2(^2) contains all the elements of IA2 which 
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contain f(x), including H^- On the other hand, since diam |.Fi([/i)| < r, every 
Wi £ Fi{U\) is contained in B(f(x) 7 r) C W%. This proves that indeed Ti is 
coarser than T\. □ 

The statement in the Theorem now follows immediately from Lemma 12.251 and 
the definitions. □ 

Remark 2.26. Lemmas 12.241 and 12.251 show a tremendous advantage of using open 
covers as opposed to partitions, because there is no easy condition defined in terms 
of the features of a partition itself which would guarantee that this partition and any 
map on it are coarser than some other partition and some other map, respectively, 
both possibly given a priori (e.g., as a result of some computation, as we do in 
Section [5]). 

Remark 2.27. If one prefers to do computations for partitions instead of using open 
covers, then J-q in the assumptions of Theorem [5] can in fact be a multivalued map 
on a partition. Definition 12.101 of the outer resolution in this case can be applied 
directly, Definition ^. 9l of a representation can be taken the same and is equivalent to 
requesting that f(U) C int |.F(Z7)|, and the predicate V can also be the same as for 
a combinatorial map on a cover. Note, however, that the conclusion of Theorem [3 
is not valid for partitions, because their inner resolution is 0. 

3. Transitivity and mixing 

Let J- : U — o U be a combinatorial map, and let J 7-1 denote its inverse defined 
by the following condition: V £ J-" _1 (C/) if and only if U £ F(V). We introduce 
the notions of transitivity and mixing for combinatorial maps, further also called 
combinatorial transitivity and combinatorial mixing, respectively. 

Definition 3.1. J 7 is transitive if for every U, V £ U there exists n > such that 

V £ J-~ n (U). J- is mixing if for every U, V £ U there exists N > such that 

V £ T~ n {U) for all n > N. 

We shall prove the following 

Proposition 3.2. Combinatorial transitivity and combinatorial mixing are finite 
resolution properties. 

To prove the proposition we need to prove the computability condition (A) and 
the consistency condition (B) in the definition of finite resolution property. Con- 
dition (A) says that the properties of (combinatorial) transitivity and mixing are 
algorithmically verifiable in finite time. This is almost immediate from the def- 
initions and Lemma 13.51 below, since combinatorial transitivity and mixing only 
require checking a finite number of conditions (because U is finite); however, in 
Section |4] we provide explicit algorithms which can be used to verify the conditions 
in an efficient way in practice. In this section we prove (B). 

Notice first that mixing implies transitivity but the converse is not true in gen- 
eral. Moreover, since V £ J 7 ~ n (U) is equivalent to U £ F n (V), one can equivalently 
formulate the combinatorial conditions for the n-th iterate of T instead of the n-th 
preimage of J 7 , thanks to the symmetry between U and V in these conditions. 
We shall use this easy observation transparently in the sequel. We shall consider 
transitivity and mixing separately. 

Lemma 3.3. Combinatorial transitivity is a consistent property. 

Proof. Let U\ and U2 be two covers of X, with U\ finer than U^- Let T\:U\ -<> 
U\ and Ti'-Ui —° M2 be combinatorial maps. Assume that T\ is finer than Ti 
and that T\ is transitive. We shall show that under these assumptions also T<z is 
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transitive. Let U2,V2 £ ^2- Let B be a ball contained in [/ 2 and disjoint from all 
the other elements of U2 ; the existence of such a ball follows from the fact that 
U2 is essential. Let U\, V\ £ be such that J7i n B 7^ and Vi C V2. Obviously, 
in this case U\ C C/2- By the assumption on transitivity of T\, there exists n > 
such that C/i <= 7?(Vi). Since C ^(V^l, we know that C/i C \F%{V 2 )\. 

In particular, PiflBc l-TO^)]; an( i therefore € ^(V^), because otherwise no 
part of B could have been covered by ^^(V^)- Since the choice of U2, V2 £ was 
arbitrary, this proves that T% is transitive. □ 

Lemma 3.4. Combinatorial mixing is a consistent property. 

Before proving Lcmma l3.4l we first prove two simple lemmas regarding an equiv- 
alent condition for mixing of combinatorial maps. 

Lemma 3.5. Let J- ': 14 —°14 be a combinatorial map, and let k > 0. If J- k (U) = 14 
for all U £14 then also !F n (U) = 14 for all U £ 14 and for all n > k. 

Proof. This follows by induction from the fact that F k+1 (U) = \J V £F(u) F k (V) 
and that F{U) ^ and F k {V) =14. □ 

Lemma 3.6. Let T: 14 — ° 14 be a combinatorial map. Then T is mixing if and 
only if there exists k > such that J rfe ([7) = 14 for all U G 14 . 

Proof. Let us first assume that J- is mixing. Then for every U,V £ 14 there exists 
Nu,v sucri that V € J- n (U) for all n > Nu,v- Since 14 is finite, the number k := 
m&x{Nu v '■ U, V € 14} is well defined and finite. Obviously, for this k, the set 
T k (U) contains every V £ 14, for all U £14. 

Now suppose that F k {U) = 14 for all U £14 and some k > 0. Lemma 1531 implies 
that J-" n (J7) = 14 for all C/ £ 14 and all n > k, which immediately implies that T is 
mixing. □ 



Proof of Lemma \3.4\ Let 14\ and U2 be two covers of X, with 14\ finer than 142- 
Let T\:14\ -<> 14\ and T2 ■ U2 —° U2 be combinatorial maps. Assume that T\ is 
finer than T2 and that T\ is mixing. We shall show that under these assumptions 
also J-2 is mixing. Take k > such that T k (U) = 14\ for all U £ 14\, given by 
Lemma [331 Take any W £ U2 and consider T\ {W). Since 14\ is finer than U2, there 
exists some U £ 14\ such that U C W. Since T\ is finer than T%, it follows that 
\T k (U)\ C \ J~2 (W) \. The assumption that U2 is essential implies that T k {W) = U 2 . 
Since the choice of W £ U2 was arbitrary, it follows from Lemma 13.61 that T2 is 
mixing. □ 

4. Graph algorithms 

In this section we provide explicit algorithms for the verification of combinatorial 
transitivity and mixing. In particular, we give a constructive proof of (A). The first 
step is to translate the properties defined in Section [3] into the language of graphs 
associated with combinatorial maps. 

We recall that a finite directed graph, further called graph for short, is a pair 
G = (V, £), where V is a finite set whose elements are called vertices, and £ c V x V 
is a set of selected (ordered) pairs of vertices. The elements of £ are called edges. 
Combinatorial maps are very naturally encoded as graphs. 

Definition 4.1. We say that G = (V, £) is the graph associated with a combina- 
torial map J" on a cover U of X if V = U and £ = {{U, V) £ U x U : V £ T{U)}. 

In this section we explain how the verification of the properties of combinatorial 
transitivity and mixing can be reduced to the verification of certain properties of the 
associated graphs. We then describe in some detail specific algorithms that verify 
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these properties, and we upper provide bounds for their effectiveness. Before we 
show such algorithms, let us explain how their effectiveness is measured. A running 
time estimate of an algorithm is typically given by means of the order of the number 
of required primitive operations as a function of the size of the input, here given 
by |V|, the size of the vertex set, and \£\, the size of the edge set. We use the 
notation 0(\&(|V|, |£|)J to indicate that there exist constants c, tiq > such that 
for any graph G — (V, £) for which |V|, \£\ > no, the number of operations $(G) of 
the algorithm applied to the graph G satisfies the inequality §{G) < c$(\V\, \£\). 
This gives an asymptotic upper bound for the worst case running time. For a more 
detailed explanation of this notation and of running time in general, the reader is 
referred to fill §3.1]. In what follows, we say that an algorithm runs in linear time 
if its worst case running time is 0(|V| + \£\). 

Proposition 4.2. There exists an algorithm which verifies in linear time whether 
a combinatorial map is transitive or not. There exists an algorithm which verifies 
in linear time whether a combinatorial map is mixing or not. 

Notice that the size |V| of the vertex set is exactly the number of elements of the 
cover U. Moreover, in many situations, such as if the map / is Lipschitz and the 
elements of U are of similar size, the expected size J-(U) of the image of each element 
U, is uniformly bounded, independently of the resolution of the cover. Therefore, 
in practice, the size \£\ of the edge set is bounded by a constant multiple of the 
size |V| of the vertex set, which is exactly the number of elements in the cover. 
Thus the algorithms we describe, in many situations, actually run in linear time in 
the number of the elements of the cover IA. We consider this linear time property 
as an important characteristic of our method since it implies that as computing 
memory and power increase, it will be realistically possible to obtain significant 
improvements to the scale at which mixing is verified in systems of interest, such as 
the Henon map which we consider in this paper (see, however, additional remarks 
on this point in Section r5.2.3[) . 

The remaining part of this section is devoted to the proof of Proposition ^. 21 In 
particular we shall obtain completely explicit forms of the required algorithms. In 
Section [4. II we show that the problem can be formulated in terms of the verification 
of certain properties of directed graphs, namely strong connectedness and aperi- 
odicity. The algorithms for verifying these properties are not new. An algorithm 
that computes strongly connected components of a directed graph in linear time 
belongs to the canons of graph algorithms (see [Til §22.5]), and we use it below to 
determine whether a given graph is strongly connected. An algorithm for proving 
aperiodicity is given in [15] but without a formal proof of its correctness, which we 
shall give here. For completeness, and for the benefit of readers who are not familiar 
with graph algorithms or algorithms in general, we shall give full descriptions of 
these algorithms below. They actually require some non-trivial constructions which 
we believe are of independent, albeit relatively technical, interest. In particular, we 
provide a fully self-contained proof of Proposition 14.21 

In Section [4.21 we discuss the representation of graphs in computer's memory. In 
Section 14.31 we describe a basic way of "exploring" graphs algorithmically and of 
describing certain structures of graphs. This approach, called "depth first search" , 
underpins both the method for computing strongly connected components described 
in Section 14.41 and the method for verifying aperiodicity described in Section 14.51 

4.1. Strongly connected and aperiodic graphs. We now show how the notions 
of combinatorial transitivity and mixing can be reduced to certain properties of the 
associated graphs. A path in a graph G = (V,£) is a sequence (Ui)f =Q such that 
(Ui,U i+ i) e £ for alH = 0, ...,n- 1. 
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Definition 4.3. We say that a graph G = (V, £) is strongly connected if for every 
U, V £ V there exists a path (vq, . in G such that vq = U and Vk = V. 

The length of a path (L^)" =0 is n. A cycle in G is a path in which Uq — U n . The 
greatest common divisor of the lengths of all the cycles in a graph G = (V,£) is 
called the period of G. 

Definition 4.4. We say that a graph G = (V, £ ) is aperiodic if its period equals 1. 

Lemma 4.5. j4 combinatorial map is transitive if and only if the associated graph 
is strongly connected. A combinatorial map is mixing if and only if the associated 
graph is strongly connected and aperiodic. 

Proof. Since the existence of a path of length n from U to V in the graph associated 
with a combinatorial map T is equivalent to V G J- n (U), the equivalence in the 
case of transitivity follows immediately from the definitions. 

The statement on mixing can be derived, with certain effort, from |17j . but for 
the sake of completeness we provide a detailed proof. Assume that a combinatorial 
map J 7 : U — ° 1A is mixing. Then Lemma 13.61 implies that there exists k > such 
that F k (U) — U for all U € U. It is thus obvious that G is strongly connected. 
Lemma 1331 implies that also F k+1 (U) = U for all U £ U. In particular, U £ T k {U) 
and U £ J- k+1 (U). In terms of the graph G, this implies that there exist cycles of 
the co-prime lengths k and k + 1 in G through U, and thus G is aperiodic. 

Let us now focus on the opposite implication. Let T = {Uj a , . . . , Uj t ) for some 
t > be a cycle in G that runs through all the vertices of G (its existence follows 
from the strong connectedness of G) . Since the GCD of the lengths of all the cycles 
in G is 1, there exist cycles C%, . . . , C r and D%, . . . ,D S of lengths p%, . . . ,p r and 
qi, . . . , q s , respectively, such that p\ + ■ ■ ■ + p r — q\ — ■ ■ ■ — q s = 1. Consider the 
cycles Tki in G composed of T, k copies of C\, . . . , C r , and I copies of D\, . . . , D s . 
The length of each such cycle is tki — t + kp + Iq, where p :— p\ + • • • + p r and 
q := q\ + ■ ■ ■ + q s . Since p = q + 1, we have tki = t + (fc + l)q + k. If we fix any 
vertex U of G at the cycle T and consider k = 0, . . . , t — 1 with I = t — k, then we 
obtain paths from U to U whose lengths are t + tq + k = tp + k. We can complement 
these cycles by following T to paths of length tp + t = t(p + 1) that end at each 
subsequent vertex of T. Since T runs through all the vertices of G, it follows that 
J rtl J>+ x ) (JJ) = U. Since U € U was chosen arbitrarily, it follows from Lemma [3.61 
that J- is mixing. □ 

4.2. Representation of graphs. The first step in developing algorithms for study- 
ing properties of a graph is to represent such a graph in computer's memory. One 
standard way to do this is as a collection of adjacency lists. We say that a vertex v 
is adjacent to another vertex u if the edge (u, v) is in the graph. The adjacency-list 
representation of a graph G — (V, £) consists of an array of |V| lists, one for each 
vertex in V. For each u £ V, the adjacency list of u contains all the vertices v £ V 
such that (u, v) £ £ in certain order. Searching through this list or taking one ad- 
jacent vertex after another can be done in the time proportional to the number of 
vertices in the list. 

Other representations, such as an "adjacency matrix" representation, can also 
be used and may be more or less convenient depending on the structure of the 
graph under investigation. As a general rule, the adjacency-list representation is 
preferable for sparse graphs, that is graphs for which \£\ is of a lower order than 
its maximum possible value of |V| 2 , e.g., if it is proportional to |V|. This is usually 
the case for combinatorial maps and thus we shall use this representation of graphs 
here. 
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4.3. The Depth First Search algorithm. Depth First Search (DFS) is an al- 
gorithm for scanning a graph G = (V,£). It starts at an arbitrarily chosen vertex, 
and follows edges and paths to visit other vertices in the graph. If no more ver- 
tices can be reached in this way then an arbitrary unvisited vertex is chosen and 
the algorithm continues running from that vertex. Each vertex u € V is assigned 
a unique discovery time t u E N which reflects the order in which the vertices are 
visited, or discovered. The general strategy of the algorithm while visiting a vertex 
is to take the first edge leading to an unvisited vertex and to follow it. If no more 
vertices can be discovered from a given vertex, then the algorithm traces back and 
checks the other edges leading from the previous vertex. The name "depth first 
search" comes from the fact that the search goes as deep into the graph as possible 
using the first edge leading to an unvisited vertex, and the other edges emanating 
from each vertex are only checked afterwards. This is different from another well 
known strategy, called "breadth first search", in which all the edges leading from 
each vertex are checked before the search continues from the discovered vertices. 
The details of the DFS procedure are summarized in the following 

Algorithm 4.6. 
function DFS 

input: G = (V, £) — directed graph 
code: time := 

for each u E V do 

if not defined t u then 
DFS-Visit (it); 

function DFS-Visit (u) 

time := time + 1; 
t u time; 

for each v E V such that (u, v) E £ do 
if not defined t v then 
DFS-Visit («); 

This algorithm runs in linear time 0(|V| + \£\), because the function DFS-Visit 
is called exactly once for each vertex in the graph, and all the edges emanating from 
each vertex are also checked exactly once. Note that in order to achieve the linear 
time in practice, it is necessary to represent the graph in such a way that scanning 
through the set of all the edges emanating from a given vertex can be done in the 
time proportional to the number of these edges, like in the case of the adjacency- list 
representation of a graph. The reader is referred to [TTJ §22.3] for more details on 
the DFS algorithm. 

4.4. Strongly connected components. One of many applications of the depth 
first search algorithm is a method for finding the strongly connected components 
of a graph. 

The most standard approach (see [TTJ §22.5]) is to run the DFS algorithm on 
the graph G, and then to compute depth first trees of G T (the transpose of G, 
obtained from G by inverting the direction of all the edges), which turn out to form 
the strongly connected components of G. Although this algorithm runs in linear 
time and memory, it is inconvenient if applied to huge graphs because of the need 
to compute G T , which takes up as much memory as G itself, and thus effectively 
doubles memory usage. 

Because of this reason, we use Tarjan's algorithm instead for computing strongly 
connected components of a graph (see [4*3114*5] ) . This algorithm is based directly upon 
DFS. In addition to computing the discovery time t u for each visited vertex, it also 
computes the lowest discovery time l u of all the vertices reachable from that vertex. 
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This number is computed during the DFS itself, so that it is always known when 
needed. Moreover, during the DFS run, the visited vertices are put on a stack S, and 
each strongly connected component Ck is taken from the stack whenever a visited 
vertex is determined to "close" such a component. For the sake of completeness, we 
provide a pseudocode of this algorithm below, but we refer to [33] for the proof of 
its correctness and for more explanations. 

Algorithm 4.7. 
function SCC 

input: G = (V, £) — directed graph; 
code: time := 0; k := 0; S := empty stack; 
for each u E V do 

if not denned t u then 
tarjan (u); 
return (Ci, . . . , C fe ) 

function tarjan (u) 

time := time + 1; t u := l u := time; 
put u on the top of the stack S; 
for each v 6 V such that (u, v) 6 £ do 
if not denned t v then 
tarjan (v); 
l u := min(7 u , l v )\ 
else if v is in S then 
l u := min(; u ,t ll ); 
if t u = l u then 
k := Jfe + 1; 

remove the elements from the top of S until u has been 

removed, too, and put all of the removed elements into Ck ; 

With this algorithm, we now determine whether a given graph is strongly con- 
nected in the most obvious way possible. Namely, we compute the strongly con- 
nected components and check if there is only one such component, as made precise 
below. 

Algorithm 4.8. 
function StronglyConnected 
input: G = (V, £) — directed graph; 
code: (Ci, . . . ,C fe ) := SCC (G); 

if k = 1 and card C\ = card V then 
return true; 

else 

return false. 

The following features of the above algorithm are obvious: 

Lemma 4.9. Algorithm \4--S\ applied to a directed graph G = returns true if 

and only if G is strongly connected. This algorithm runs in linear time. 

4.5. Aperiodic graphs. Efficient computation of the greatest common denomina- 
tor of cycles in a strongly connected directed graph is less standard, and the solution 
suggested in [17] yields cubic time. However, there exists an algorithm which runs 
in linear time, see [JS]. For the sake of completeness, we describe this algorithm 
below and prove its correctness (note that the latter is not provided in [48 ). This 
algorithm is based upon the standard DFS (Depth First Search) algorithm that 
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scans the entire graph starting at an arbitrarily chosen vertex (recall the assump- 
tion that the graph is strongly connected), and computes the greatest common 
divisor of certain numbers related to the edges of the graph. 

Algorithm 4.10. 
function GraphPeriod 

input: G — (V,£ ) — a strongly connected directed graph; 
u — an element of V; 

d u G Z — the depth of u in the DFS tree being constructed; 
p G Z — the GCD of cycle periods found so far; 
code: for each v G V such that (u, v) G £ do 
if defined d v then 

p:= GCD (p, d v -d u -l); 

else 

GraphPeriod (G, V, d v := d u + 1, p); 

return p. 

Lemma 4.11. Algorithm ^. 1Q\ applied to a strongly connected directed graph G — 
(V, £), any element r G V, and the numbers d r := and p := 0, returns the greatest 
common divisor of the lengths of all the cycles in G. This algorithm runs in linear 
time. 

Proof. Given a strongly connected graph G = (V, £), let T = (V, £') denote the tree 
obtained by running Algorithm 14.101 on G, starting with the vertex r G V, which 
becomes the root of T. The set £' consists of all the edges which incur recursive 
calls of the function GraphPeriod. The computed depth of each vertex u G V in T is 
denoted by d u , with d r = 0. For each edge e = (u, v) G £ \£' , define d e := d v — d u — l. 
Let c be the GCD (greatest common divisor) of all these numbers d e . This is the 
quantity returned by this algorithm. We shall prove that c equals the GCD of the 
lengths of all the cycles in G, further denoted by d . 

Let us first prove that c'\c. Consider e = (u, v) G 8 \ £' . Since G is strongly 
connected, there exists a path p v>r in G from v to r; denote its length by d. Note 
that the path p rA1 in T from r to u is of length d u , and the path p r ^ v in T from r 
to v is of length d v . Since p r>u combined with e and p v<r is a cycle in G of length 
d u + 1 + d, and p r>v combined with p v ^ r is a cycle in G of length d v + d, the common 
divisor d of the lengths of all the cycles in G must divide the difference between 
these lengths, which is d v — d u — 1. Therefore, c is a GCD of numbers divisible by 
d, and thus d\c. 

Let us now prove that c\d . Let (uq, . . . , v n ) be any cycle in G. We shall prove 
that c\n. Consider the depths of each Vi in T. The difference in the depth traversed 
by each edge = (vi, i>j+i) is Si := d Vi+1 — d Vi , which equals 1 if e; G £' . Since wo = 
v n , obviously d Vo = d Vn , and thus Y^i=o ^ = 0- ^ n particular, — n = Yll=o 1) = 
S™=o 1 (^"i+i — ^» — ^ e i ^ ^' t nen f ne corresponding item in the sum is zero, 
otherwise it is divisible by c. Therefore, c\n. 

The observation that each edge in the graph G is processed exactly once shows 
that the algorithm runs in linear time 0(|£|), which completes the proof. □ 

5. Numerical computations 

We are finally ready to apply all the ideas above to a specific example. First of 
all, we state a more precise version of the theorem given in the introduction. 

Theorem 3. There exists an open set X C IR 2 such that H a ^(X) C X and H a ^\x 
is (combinatorially) mixing at all resolutions > 1CP 5 for all (a, 6) in an open set 
Pel 2 containing (1.4,0.3). 
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We choose the Henon map for the so-called "classical" parameter values simply 
because it is a very well known and well studied example for which nothing sub- 
stantial is known. We have purposefully formulated all our definitions and results so 
far in full generality so that it is straightforward how to apply them to many other 
examples. What is left is to construct a representation T of the Henon map such 
that TZ + (J-) < e and then apply the algorithms described above to show that T 
is mixing. We describe the construction of J- in detail so that minor modifications 
should allow essentially the same construction to work in many other cases. 

Notice that the result applies to an open set of parameter values. This is an in- 
trinsic feature of the method deriving from the use of interval arithmetic. Since both 
numbers 1.4 and 0.3 are not representable in the binary floating-point arithmetic, 
small intervals containing these numbers are taken for the actual computations. 
Therefore, the statement is automatically proved for an open set of parameters 
containing these numbers. This costs nothing from the computational point of view 
and in fact corresponds to the natural fact that numerical methods necessarily only 
yield results which are stable under sufficiently small perturbations. 

5.1. General strategy. We describe in detail the strategy for a general continuous 
map /: R™ -> R". 

5.1.1. Open interval arithmetic. For the purpose of the numerical computations, 
we use a slight modification of interval analysis introduced in |29j . 

In the typical approach, closed intervals are used instead of numbers, and the 
result of each arithmetic operation on such intervals is defined as the smallest possi- 
ble interval containing the results of the operation on any numbers taken from each 
of the intervals, e.g., [01,02] — [61,^2] = [0*1 — 02,0-2 — &i]- Since we are interested 
in doing these computations using a fixed-size floating point representation of real 
numbers, the set of actual numbers that can be represented is finite, and thus we 
must round the endpoints of the resulting intervals to the nearest representable 
number in the downwards direction for the left endpoint, and in the upwards direc- 
tion for the right one. As a consequence, the result of calculations carried out on 
intervals is an interval that contains every possible exact result of those operations 
on the numbers belonging to these intervals. 

Since we work with open covers, in our case it is necessary to work with open 
intervals, instead of closed ones. This implies slight differences in the case of some 
arithmetic operations on intervals. Although addition, subtraction, multiplication 
and division are the same, those operations in which closed intervals would arise 
must be slightly changed, like rising to even powers (e.g., computing x 2 ) or some 
trigonometric functions (mainly sin and cos). Namely, the result of an operation on 
open intervals with representable endpoints is defined as the smallest open interval 
with representable endpoints which contains the results of the operation on single 
elements taken from these intervals, e.g., (—1, l) 2 = (— e, 1), where — e < is the 
largest representable negative number (the set of representable numbers is finite, 
so this number is well defined). 

5.1.2. Open cover parameters. We start by selecting a bounded rectangular area 
B := (ax,ai + wi) x ■ ■ • x (a„,a„ + w n ) C R" which is assumed to contain the 
dynamics of our interest. We cover this area with a finite family of overlapping 
open boxes (n-dimensional hypercubes). We set some number p\ £ N of parts into 
which we intend to subdivide the first interval (ai,a\ + W\) in the definition of 
B, and then we compute the numbers p%, . . . ,p n so that each pi is approximately 
proportional to Wi for i = 1, . . . , n. Then we choose some small n > to be the size 
of the overlapping margin. As a general rule we can choose k to be the smallest 
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positive representable number in the chosen computer precision (recall that the set 
of representable numbers in (0, 1) is finite). We then consider the family 

B :« {n™ =1 (ai + hwi/pi,^ + (ki + l)wi/pi + n) : ki = 0, . . . ,Pi - l}, 

where the :« symbol is used to indicate the fact that the actual endpoints of the 
intervals are computed in the floating point arithmetic with rounding the result of 
each operation to the nearest representable number either in the downwards direc- 
tion (for the left endpoints) or in the upwards direction (for the right endpoints). 

5.1.3. Construction of X. The final cover hi, and thus the actual space X, is chosen 
by an iterative procedure from the elements of B, so that U C B. In particular, this 
cover is necessarily essential. 

The iterative procedure is carried out as follows. We start with a point xq € 
K™ as close to the attractor as possible. This point may be obtained for example 
by some non-rigorous numerical simulation of the dynamics, but the construction 
is in general not particularly sensitive to this choice. We then define an initial 
approximation of the cover by letting 

Ua := {U e B : x e U}. 

For each (7 £ Wo we compute an open set f(U) as the image of U under / using 
interval arithmetic. This is a rigorous upper bound for f(U) in the sense that 

/(CO c W). 

We then define the multivalued map J-q : Uq — ° B by 

T Q {U) :={BEB:BDf(U)^(t)}. 

At this point we compare Uq with Tq(Uo) = \JueU -Fo(U). If J 7 o(^o) C Uq then we 
let U := Uq and define J 7 := -Fo|w : U -^U. This is our combinatorial representation 
of the map /: X — > X on the set X := \U\. If Fq(Uq) <£_ Uq then we define 

and we compute the map T\\U\ —° B as an extension of J-'q using interval arith- 
metic. Note that only the images of the sets in J r o(Wo)\Wo have to be computed. Pro- 
ceeding again in the same way, we repeat this procedure and we obtain Uq, . . . , Uk, 
as well a,s J-q, ... ,J-k, until we eventually get J-(Uk) cUk, at which point we define 

U:=U k and T := T k \ u : U -<> U. 

Although this process is guaranteed to terminate, because of the finiteness of B, 
the obtained result might be faulty if /(C) (£ B for some U 6 W, as then T is 
not a combinatorial representation of f\x- If this happens then we say that the 
construction fails, and one has to choose different parameters for the open cover 
(Section I5.1.2[) and try again. In particular, if the system has an attractor with 
a sufficiently large basin of attraction and strong enough convergence, then the 
computations should result in a valid map T provided the set B is taken large 
enough, the number pi is high enough, and the accuracy of computations is good 
enough. 

Note that if this construction succeeds then the open set X :— \U\ is contained in 
the attraction basin of the attractor, f(X) c |.F(A)| c X, and J 7 is a combinatorial 
representation of f\x '■ X — > X. 
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5.1.4. Verification of the properties of J 7 . The outer resolution TZ + (J 7 ) of the map T 
can be easily calculated using interval arithmetic during the computation of J-'i, 
based on the diameter of a cover of f(U) whenever Ti{U) is constructed. 

After the combinatorial map JF has been constructed, one can verify whether it is 
transitive and mixing by a straightforward application of Algorithms 14.81 and 14.101 

5.2. Application to the Henon map. We now give the specific data involved in 
the calculations and discuss in more details the numerical and computational issues 
involved. 

5.2.1. Constructing a combinatorial representation. We closely follow the general 
strategy described in Section I5~T1 We choose the rectangular region B in such a way 
that an approximation of the attractor found in numerical simulations is contained 
in B with some safety margin (about 0.1 at each side). More precisely, we set 



where the symbol indicates that we take a representable number close to the 
given decimal number (it is most likely the nearest representable number, depending 
on a particular compiler, but it does not really matter). Note that the numbers 
which are powers of 2 are representable and thus the actual values of a 2 and u>2 
used in the calculations are exact. We choose k to be the smallest positive (normal) 
number representable in the standard double precision floating point arithmetic, 
which is approximately 2 • 10 -308 . 

The number p\ is a parameter of the program that can be changed at each run 
(it is set up from the command line each time the program is launched), and we 
figured out by trial and error that for p\ := 2,132,419 (for which the corresponding 
Pi was taken as 761,578) the algorithm constructs a map J 7 with 1Z + (J 7 ) < 10 -5 . 

As an initial point supposedly very close to the attractor we take 



This point was found in numerical simulations by iterating the origin a little over 
100 million times. 

The parameters of the Henon map for which we want to carry out the com- 
putations, a = 1.4, b = 0.3, are not representable in the binary floating point 
arithmetic. Therefore, we take open intervals containing them, computed in the 
program as the quotients 14/10 and 3/10, respectively, with the left endpoints of 
the intervals rounded down to the nearest representable numbers, and the right 
endpoints rounded up. 

For each U € B, which is a product of two open intervals, we calculate f(U) by 
means of open interval arithmetic, as described in Section [5 .1.11 using the formula 
for the Henon map and the intervals containing the parameters a and 6, as explained 
above. 

Instead of constructing subsequent maps J-'i, the algorithm is set up in such a 
way that it constructs a list U containing the elements of the initial cover Uq (which 
covers the initial point xq), and then for each element of this list it computes F(U) 
and immediately appends the elements of FQJ) which are not yet in the list U to 
the end of this list. The algorithm terminates if the end of the list has been reached, 
which means that all the elements of F(U) for every U in the list are already in 
the list. The algorithm quits with a failure result if f(U) <f. B. A sample result of 
the constructed set U is illustrated in Figure 15.2.11 



ai :« -1.4, 
a 2 := -0.5, 



wi :« 2.8, 
w 2 := 1.0, 



x :« (0.61989426930989,0.17586130934794). 
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Figure 3 . A sample low- resolution cover of the Henon attractor 
with close-ups of selected details. 

5.2.2. Applying the graph algorithms. The constructed cover IA and combinatorial 
map T are represented in two data structures: an array of pairs of open intervals 
which represent the elements of U (that are products of these pairs of intervals), 
and a directed graph whose vertices are the indices to this array. After the map 
JF has been constructed, the memory used for the array representing hi is released, 
and the graph algorithms are run on the graph representation of J- ' . 

In the actual implementation of Algorithms 14.71 and 14.101 instead of using the 
recursive call to the subroutines DFS-Visit and tar j an, respectively, we use a 
stack version of DFS because of the limitations of some systems which allow for 
a limited recursion depth only. In this version, instead of calling the subroutine, 
the current parameters are put on the stack and another round of the main loop 
is started with new values, with a return from the recursive call corresponding 
to taking the previously stored values of the variables from the stack. This is a 
standard technique, and one can check the details of its implementation in our 
software available at [55] . 

5.2.3. The cost of the computations. As it should be expected, the time and memory 
usage of the computations for the Henon map heavily depend on the number p\ . For 
small values of pi, up to some 20,000, the time of computation on a computer with 
a contemporary modern processor (Intel® Xeon® 5030 2.66 GHz was used in our 
computations) should not exceed 10 seconds and use a negligible amount of memory 
(up to about 25 MB). The smallest number p\ for which the computations were 
successful was pf m = 446 (with the corresponding pf 111 = 159, and the computed 
7Z + (J- mm ) was below 0.05. For smaller values of p\ (and also for a few larger ones) 
the algorithm fails because of encountering the situation in which f(U) <f_ B for 
some U in the cover being constructed. 

For higher values of p%, we observed that if p\ is increased 10 times then the outer 
resolution of the computed map T decreases about 10 times, while the number 
of elements in the constructed cover U grows about 20 times, and so do both 
the computation time and memory usage (see the actual results available at (3"5]b 
In particular, for the highest tested p™ ax ;= 2,132,419 (with the corresponding 
^max _ 761^578), a cover £/ max consisting of 161,448,094 boxes was constructed, 
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and 7?. + (J rmax ) was slightly below 10 -5 . The computation time was about 1 hour 
and 17 minutes, out of which the majority was used for the construction of the cover 
U and the graph associated to the map T , while running the graph algorithms took 
less than 6 minutes. 

In these computations we reached the limits of the computer equipment available 
to us at the time of writing this paper, but with the development of technology 
one can speculate on the possibilities of doing more extensive computations in 
the future. Based on our rough estimates, obtaining the resolution 10~ 6 would 
require some 26 hours of processor time and about 360 GB of memory (twice 
more than 20 • 9 = 180 GB because of the need to switch from 32-bit to 64- 
bit integers in the graph representation, due to the high number of vertices and 
edges). Getting down to 10 -7 would take over 500 hours (3 weeks) and use some 
7 TB of memory. Without switching to distributed computations and probably also 
changing the approach (e.g., re-computing the graph on-the-fiy instead of storing it 
in the memory) , these and higher resolutions seem to be still out of reach for many 
years to come. However, one might argue that such high precision of the results is 
not necessary for real applications, and 10 -5 is more than one might require, so 
investing in better results does not make sense at this point. 

Finally, we would like to point out that the method we have developed is dimension- 
independent, and so is the related software available at |38| . However, the number 
of elements in a typical cover of a complicated attractor should be expected to grow 
considerably faster in higher dimensions with the increase of the target resolution. 
This might be a major limitation of the applicability of our method for obtaining 
results at very fine resolutions for more complicated systems. In practice, however, 
even complicated dynamics in high-dimensional spaces often concentrates around 
low-dimensional attractors (like the topological horseshoes), in which case the res- 
olutions interesting from the point of view of applications might be still within the 
reach of contemporary computers, and we sincerely hope to see such applications 
in the future. 

6. Final remarks: Finite resolution versus topological and 
measurable dynamics 

It is tempting to try to see finite resolution properties as approximations of 
the "real", for example topological, properties of the underlying map f:X—>X. 
In terms of the transitivity and mixing properties considered here, this is partly 
justified in the sense that if / is topologically transitive or mixing then it is also 
combinatorially transitive or mixing, respectively, at all resolutions (for complete- 
ness, we give a formal proof in Appendix |B|) . However our combinatorial notions 
of transitivity and mixing are not uniquely defined by this requirement, and to 
fully justify the definitions and the idea that we are approximating the topological 
properties we would need to show the converse result that if / is combinatorially 
transitive or mixing at all resolutions then it is topologically transitive or mixing. 
It is possible to obtain such a double implication for certain kinds of properties, see 
for instance the comprehensive discussion on precisely this topic in |30j . However, 
as mentioned in the introduction, this can be achieved only for "robust" properties 
which are persistent under small perturbations of the system. It seems therefore 
that in general the finite resolution dynamical properties of a system cannot be 
thought of as an approximation of the topological properties, at least not in a naive 
sense. 

For a time we struggled with this limitation of the theory and perceived it as a 
weakness. However, on further reflection we realized that it was in fact arising from 
an incorrect understanding of what this approach is actually about. The mistake 
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was to assume that the topological point of view is in fact the "real" dynamics and 
thus the ultimate goal. In fact, on a fundamental level, the real dynamics, if we really 
want to have such a notion, has to be just the dynamics of / thought of as a function 
on the set X. Any more sophisticated description necessarily relies on additional 
structure and the description then necessarily has to be in a form which is com- 
patible with this structure. There are at least two major frameworks or structures 
through which, in many cases, one and the same dynamical system can be studied: 
the topological and the measurable, which provide in some sense two alternative 
points of view on the same dynamics. Each of these points of view comes with its 
own definitions, notions, and results which are intrinsically motivated within that 
particular framework. A good example, very relevant to the present paper, is pro- 
vided by the notions of topological mixing and measure-theoretic mixing. Both of 
these notions formalise some intuitive notion of "mixing" within the corresponding 
framework, but there is no direct relationship or formal implication between them 
in general (some systems may be topologically mixing but not measure-theoretically 
mixing and vice versa) . Analogously, the notion of mixing which we give in the pa- 
per is a combinatorial notion which is very natural in the finite resolution setting. 
The relationship of this notion with those of topological and measure-theoretical 
mixing is interesting and quite complex, and certainly deserving of further thought, 
though beyond the scope of this paper. 

In conclusion, we believe that the framework of finite resolution dynamics pro- 
vides one of several possible structures through which to study the dynamics. Ac- 
cordingly, the definitions and dynamical features of interest should be intrinsically 
motivated within the finite resolution framework. Thus finite resolution dynamics 
should be seen as an alternative structure, alongside the topological and the mea- 
surable, which can contribute to an effective study of a dynamical system from a 
different point of view. 

Appendix A. Essential covers 

Proposition A.l. From any cover it is possible to create an essential one (without 
increasing 1Z + , but possibly decreasing the number of elements). 

Proof. Let U° = ([/{*, . . . , be an ordered cover that is not necessarily essential. 
We shall construct a series of gradually "corrected" ordered covers U 1 , . . . ,U n such 
that each of the covers has exactly n elements, U k = (U k ,...,U k ), and each of 
the first k elements of U k either is empty, or satisfies the exclusive ball condition 
that appears in Definition 12.51 of an essential cover. Then U n will give rise to an 
essential cover after having removed the empty sets. We proceed by induction. Note 
that U° satisfies the inductive assumption. Suppose that U k has been constructed 
for some k < n. Then each U k for i < k is either empty or satisfies the exclusive 
ball condition, and then we set U k+1 :— U k for i < k. Consider U k +1 . If it is 
empty or satisfies the exclusive ball condition then take U k+1 :— U k also for all 
i = k + 1, . . . , n. If Ug +1 C U^fc+i u i then set u kti : = and U i +1 '■= U i for a11 
i = k + 2, . . . , n. Otherwise, take any x € U k +1 . Since U% +1 is open, there exists 
r > such that B(x,r) C U k +1 . Take U%+{ := U k +1 and U k+1 := U k \ cl B{x, r/2) 
for i = k + 2, . . . , n. Note that all U k+1 , i — 1, . . . , n, are open sets (although some 
might be empty), and U™=i U k = U"=i so U k+1 is an ordered cover of the 

same space. Since E7- 1 = U k for i < k and U k+1 C U k for i > k, the exclusive ball 
condition for each nonempty U k+1 with i < k follows directly from the one for the 
corresponding U k . Moreover, by the construction, is either empty, or satisfies 

the exclusive ball condition with B(x,r/2). □ 
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Appendix B. Topological mixing implies finite resolution mixing 

Proposition B.l. Suppose X is a metric space and f : X — > X a continuous map 
which is topologically transitive or topologically mixing. Then any combinatorial 
representation T of f is transitive or mixing, respectively. 

Proof. Recall that / / is transitive if for every two open sets U, V C X there exists 
n > such that f~ n {U)f\V 7^ 0. Let T : U — o U be a combinatorial representation of 
a transitive map / : X — > X. We shall show that T is transitive. Take any U,V eW. 
Since / is transitive, there exists n > such that f~ n (U) n V 7^ 0. Therefore, 
/"(/-"({/)) n f n (V) ^ 0. Since /"(/""(f/)) C E/, we have [/ n /"(F) 7^ 0. Note 
that T n is a representation of /", and thus U G F n {V). Since the choice of U and 
V was arbitrary, this implies that T is transitive. 

Similarly, / is mixing if for every two open sets U, V C X there exists N > 
such that f~ n (U) n V 7^ for all n > TV. Let J": W -o W be a combinatorial 
representation of a mixing map /: X — > X. We shall show that J 7 is mixing. Let 
Since / is mixing, there exists iV > such that / n (U) H V 7^ for all 
n > iV. Fix this n. Note that f n {f-' n {U)) n /"(V) 7^ 0. Since f n {f~ n {U)) C [/, we 
have £/Tl/ n (V) 7^ 0. Recall that J 7 ™ is a representation of /", and thus [/ e F n {V). 
Since the choice of U and 1^ was arbitrary, this proves that T is mixing. □ 
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